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... Abstract 
X , 

V" Heavy neutrino with rrih < 300MeV/c^ can be effectively searclied for in kaon 

^ decays. We put upper limits on a mixing matrix element |f//^/iP for radiatively decaying Uh 

^ ! from K~ — )■ ^~Vh{i-'h ^l) decay chain in the following parameter region: 30MeV/c^ < 
nih < 80MeV/c2; IQ-i^s < Th < IQ-^s. For the whole region \U^h? < 5 ■ 10"^ for 

I Majorana type of Vh and < 8 ■ 10"^ for the Dirac case. 

O ■ Keywords: radiative kaon decays, heavy neutrino, sterile neutrino, LSND anomaly 
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1. Introduction 



For more than ten years latest results of short-baseline neutrino experiments have been 
widely discussed and still there is no clear understanding of an event excess observed by 
LSND \i] and MiniBooNE [sj experiments and their contradiction with KARMEN fl] 
^ I results. 

Oscillation interpretations of the event excess require additional sterile neutrino(s) 
with Am^ ~ leVVc^ (see p] for review). An alternative interpretation of the results of 
all three experiments is proposed in [6] . Below we briefly discuss results obtained in that 
paper. 

The main idea (proposed for the first time in 0]) is that in the experiments mentioned 
above signals from electrons and photons are indistinguishable. One could introduce 
heavy sterile neutrino Vh as a component of z/^ flavor eigenstate with a corresponding 
mixing matrix element f/^^ which is produced in z/^ neutral current (NC) interactions and 
decays radiatively into a photon and a light neutrino v. The decay channel — )■ is 
dominant if there is a large enough magnetic transition moment ixtr (it requires substantial 
new physics because in a minimally extended SM is not large enough, see [sj). In this 
case the event excess in LSND and MiniBooNE experiments comes from photons and not 
from electrons. In KARMEN experiment, z//i's with m > 40 MeV/c^ cannot be produced 
within the detector because of a kinematic threshold effect. 
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Sterile neutrino could be either of a Dirac or Majorana type. In the latter case a 
photon angular distribution in rest frame is isotropic while for the Dirac case there is 
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an anisotropy depending on Uh mass: ^^g^* ~ (1 + -;^^^-^~^cos6*). 

The combined analysis of LSND, KARMEN and MiniBooNE data results in the fol- 
lowing properties of Uh (regardless of the neutrino type): 



• 40MeV/c2 <mh< SOMeV/c^; 

• 10-"s <Th< lO-^s; 
. 10-3 < |f/^,|2 < 10-2. 

It was mentioned that z//j could be considered as a component of i/^. This leads to a 
very important consequence that Uh is also produced in charged current (CC) interactions 
and can be effectively searched for in pion and kaon decays (this idea was proposed in 
i). 

The simplest way to do it is to study two-body decays tt — )■ /iz/ and K ^ fiu and look 
for a peak in the muon energy distribution (i?^ = (M^ + — m\)/2M) below the main 
one from tt — )■ /ii/^ (7r^2) and K — )■ {.K^2)- These decays allow to search for Vh with 
masses up to ~300 MeV/c^. 

Experimental limits from 7r^2 decay [loj were obtained for SMeV/c^ < ruh < SOMeV/c^: 
|f/^/jp < 10-^ — 10-3. Best limits for kaon decays come from KEK experiment ll| : 
\U^h\^ < for TOMeV/c^ < nih < SOOMeV/c^. 11^2 decay is not sensitive to large rrih 
masses, while K^2 decay is not sensitive to low masses because of resolution effects 
and strong background from K ixp^^[K^2'^) decay. Thus, a region SOMeV/c^ < rrih < 
TOMeV/c^ is not constrained at all. One should notice that limits above were obtained 
for relatively long-lived neutrinos flying away from a detector (photon veto was applied 
in both cases). 

Another possibility to search for heavy neutrino in kaon decays (which we are going 
to use) is to measure K — fiUhi^h — ^ decay chain. In this case the background from 
K^2 is small and one can search for in a low mass region (the background from 
is also small, see Section One should stress here that only the case of radiatively 
decayed neutrinos is considered. 

The main purpose of this Letter is to search for heavy neutrino in K~ — )■ ^~Vh{vh 
wy) with the properties described in j6| and in the following parameter range: 30MeV/ c^ < 
rrih < 80MeV/c2, lO-^^s < th < I0~^s. 

2. ISTRA+ setup 

2.1. Experimental setup 



The experiment was performed at the IHEP 70 GeV proton synchrotron U-70. The ex 



perimental setup ISTRA+ (fig. H]) was described in details in [l2|. The setup was located 



in the negative unseparated secondary beam. The beam momentum in the measurements 
was ~ 26 GeV/c with Ap/p ~ 1.5%. The fraction of K~ in the beam was ~ 3%. The 
beam intensity was ~ 3-10^ per 1.9 s U-70 spill. The track of a beam particle deflected by 
the magnet Mi was measured by BPC1—BPC4 (1mm step multiwire chambers), the kaon 
identification was done by Co — C2 threshold Cherenkov counters. A 9 meter long vacuum 
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Figure 1: Elevation view of the ISTRA+ detector. 

decay volume was surrounded by the Guard System (GS) - 8 lead glass rings LGi — LGs 
used to veto low energy photons. SP2 was a lead glass calorimeter to detect /veto large 
angle photons. Tracks of decay products deflected in M2 with ITm field integral were 
measured by PCi — PC3 (2mm step proportional chambers); DCi — DC^ (1cm cell drift 
chambers) and finally by 2cm diameter drift tubes DTi — DT4. Wide aperture threshold 
Cherenkov counters C3, C4 were filled with He and were not used in the measurements. 
Nevertheless C3 was used as an extension of the decay volume. SPi (ECAL) was a 576- 
cell lead glass calorimeter, followed by HC (HCAL) - a scintillator-iron sampling hadron 
calorimeter. HC was subdivided into 7 longitudinal sections 7x7 cells each. MH was 
a 11x11 cell scintillating hodoscope used to improve the time resolution of the tracking 
system, MuH was a 7x7 cell muon hodoscope. 

The trigger was provided by 5*1 — 6*3, S'5 scintillation counters, Cq — C2 Cherenkov 
counters, the analog sum of amplitudes from the last dinodes of the 5* Pi : Tq = 5*1 ■ 6*2 ■ 

■ Cq ■ Ci ■ C2 ■ ■ S(S'Pi), here S'5 was a counter downstream the setup at the beam 
focus; S(S'Pi) - a requirement for the analog sum of ECAL amplitudes to be above ~3 
GeV. The last requirement served to suppress the K^2 decay. About 10% events were 
recorded with a different trigger: Ti = Si ■ S2 ■ S-^ ■ Cq ■ Ci ■ C2 ■ S^. This prescaled trigger 
allowed to calculate the trigger efficiency as a function of the energy released in ECAL. 

2.2. Data and MC samples 

We use high-statistics data collected in Winter 2001 run. About 332M events were 
stored on tapes. This statistics was complemented by 200M MC events generated with 



Geant3 [13[. The MC generation includes a realistic description of all ISTRA-t- detectors. 

For the signal simulation kaons are forced to decay into a muon and Vh- A new particle 
Vh is introduced to Geant3: 

• Vh mass is 30, 40, 50, 60, 70 and 80 MeV/c^; 

• Vh lifetime is 10"^ 10"^° and lO'^s; 

• Vh decays into a photon and a massless neutrino; 

• photon angular distribution in Vh rest frame is isotropic (Majorana case). An 
anisotropy for the Dirac type of Vh is obtained by weighting MC events. 

Each sample of the signal MC (for particular values of uih and Th) contains IM event. 

3. Event reconstruction 

3.L Photon momentum reconstruction 

Event reconstruction for K~ — )■ ^^Vh{vh — )■ z/7) is nearly the same as for K~ — )■ 
The only difference is that for K" — it is possible to reconstruct the 



gamma momentum in the laboratory frame using the decay vertex and a shower 
centre while for — ?► — )■ 2^7) the photon is emitted from a secondary vertex 
which is unknown. Nevertheless one can reconstruct •p^"'^ using a primary vertex. This 
leads to an additional photon energy smearing in the kaon rest frame. 

In Fig. |2] the normalized difference between the measured and true photon energy in 
the kaon rest frame is shown. For r/j = 10~^^s this difference is dominated by resolution 
effects while for large Th one can see an additional smearing. Typical values of the photon 
transverse momentum with respect to Vh momentum are small, that is why the smearing 
is not crucial. 
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Figure 2: Photon energy smearing in the kaon rest frame; dE/E = {Emeasured — Etme) / Etme (signal 
MC, mft=60MeV/c2). 



3.2. Signal signatures for K — )■ iJiUhiyh ^l) 

The simplest way to observe heavy neutrino is to look for a peak in - muon energy 
in the kaon rest frame (see Section [1]). 

Another signature comes from a photon angular distribution in the kaon rest frame. 
For this distribution using Lorentz boost transformation one can obtain the following 
formula: 



dN 



dN 



1 



dcosO dcosO* -f^{l3cose - 1)^ 



where 6 is the angle between ph and p^ in the kaon rest frame, 6* is the angle between 

The term -r^^ 
dcosO* 



p* in Vh rest frame and the boost axis (along p/^), 7 



is constant 



dN 



m,, 
1 + ^ 



m 



ml + ml 



for the Majorana type of Vh. In the Dirac case ^^gg* 

In all cases 7 > 3, therefore /3 ~ 1 and the angular distribution has a peak at cos9 ~ 1 
and hence cos6^^ ~ — 1 in the kaon rest frame. This peak is a very good signature for 
the signal. 

Distributions over cosO^^ and y for Dirac and Majorana cases are shown in Fig. [31 
It can be seen that the difference in y and cosO^^ shapes is negligible for two neutrino 
types. 



3.3. Primary and secondary decay vertex 

The difference between ^-coordinates (z-axis is collinear with the beam direction) of 
the secondary and primary vertices divided by the distance between the primary vertex 
and the electromagnetic calorimeter is shown in Fig. |H A distribution over this ratio 
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Figure 3: Dirac (solid) and Majorana (dashed) type of I'h {irih—GOMcV / , Th ~ IQ^^^s). Left: distribu- 
tion over normalized muon energy y = 2Efi/mK (lines - selected region, see Section l5.2p . Right: cos9^j_^ 
in the kaon rest frame. 



shows a fraction of decays within the decay volume. For = lO^^^s and = 
lO^^'^s almost all neutrinos decay within an experimental setup, while for = lO^^s a 
geometrical inefficiency becomes substantial. 
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Figure 4: The difference between coordinates of the secondary {zh) and primary {zk) vertices divided 
by the distance from zk to ECAL (MC, m/i^GOMeV/c^). 



4. Event selection 

The event selection for — t- ^^Vhi^h ^l) is very similar to that of — )■ 
Standard kinematic variables are used for the further analysis: x = 2E^/mK and ji 



2E^/mK, and E^ being photon and muon energies in the kaon rest frame. As in 19 
Dalitz-plot will be used for studying signal and background kinematic regions. 

The decay signature is defined as follows: one primary track (kaon); one negatively 
charged secondary track identified as muon; one shower in ECAL not associated with the 
charged track. Muon identification using ECAL and HCAL is described in our previous 
papers (Qllij)- 

Several cuts are applied to clean the data: 

• number of beam and decay track projections in XZ and Y Z planes is equal to 1; 

• CL (confidence level of a track fit) for the beam track projections in both planes 
must be greater than 10~^; 
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CL for the decay track projections is greater than 0.1 (XZ) and 0.15 (YZ); 



• the angle between a primary (kaon) and secondary (muon) track is greater than 
2 mrad. 

The last cut eliminates most of undecayed beam particles. The quality of the decay track 
(described quantitatively by CL) is worse than that of the beam track because of multiple 
scattering and detector resolution. 

Cuts containing a photon energy include: 

• gamma energy in the kaon rest frame is greater than 10 MeV; 

• no photons in SP2 calorimeter (the energy threshold is 0.5 GeV for the total energy 
release) ; 

• no photons in GS. 

For vertex characteristics we have the following requirements: 

• 2;-coordinate must be within the interval 400 < z^tx < 1600cm; 

• (-3) < Xytx < 3cm; 

• (-2) < y^tx < 6cm; 

• CL of general vertex fit is greater than 10~^. 
Additional cuts are applied to suppress backgrounds: 

• number of hits in the matrix hodoscope [MH) is less than 3; 

• missing momentum "^^^^s = — — ^^^^ does not point to the ECAL central 
hole (this cut effectively rejects the background from K~ — > tt^tt^ decay when the 
lost photon from vr^ — )■ 77 goes into the hole). 

4.1. Trigger efficiency 

As To trigger described in Section 2 contains the energy threshold in S'Pi, the trigger 
efficiency as a function of energy released in ECAL should be known. It could be found 
using events with Ti trigger: Strg = (Ti f]To) / Ti. The trigger curve is shown in Fig. [5l 
The fit is done using a Fermi function. For the further analysis only events with Tq are 
kept and these events are weighted by the factor of l/strg- 

5. Signal extraction 

As it was mentioned in Section HI Dalitz-plot is used for signal and backgrounds 
studies. The main background comes from 3 decay modes: — t- /i^z/^7(_K'^27), K~ — > 
/i^z/^7r°(ir^j3) with one gamma lost from 7r° — 77 and — t- 'k^'k^{Kt^2) with one gamma 
lost and vr misidentified as /x. Dahtz-plot distributions for the signal, K^2'-i^ and Kj,2 
are shown in Figs. [6] - [91 
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Eecal, GeV 

Figure 5: Tq trigger efficiency. Points - data, curve ~ fit by the Fermi function. 




Figure 6: Dalitz-plot density for the signal Figure 7: Dahtz-plot density for the K^2-f{IB) 

(m/i=60MeV/c^, th — 10^^°s). background. 

5.1. Signal extraction procedure 

The procedure starts with dividing all the kinematic {y, x) region into stripes on x 
(x-stripes). The x-stripe width is Ax=0.05 [AE^ ~ 12MeV). In every x-stripe we put a 
cut on y: 1.0 < y < 1.2 to suppress backgrounds. 

In Section [22] two signal signatures were described - peaks in y and cos9^^. For each 
x-stripe we do a simultaneous fit of two histograms - cosO^^ (with the cut on y introduced 
above) and y (without this cut). 

5.2. Selected kinematic region 

For the further analysis we have selected seven x-stripes in the following region: 0.2 < 
X < 0.55 (49MeV < < 136MeV). Selected x-stripes are shown in Fig. |6l Dalitz-plot 
for the data with selected x-stripes is shown in Fig. [TOl 

5.3. Possible signature for different nih and tu 

To better understand how the signal looks like and how sensitive our data are to 
heavy neutrino we plot y and cosO^^ distributions for ruh = 60MeV/c^, th = 10~^°s 
and \Uf^h\'^ = 10~^ in three x-stripes (Figs. [TTl [T^ . It can be seen that cos9^^ is much 
better for the signal observation (or setting upper limits) while y is good for the reliable 
background normalization. 
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Figure 8: Dalitz-plot density for the K^-^ back- Figure 9: Dalitz-plot density for the back- 

ground, ground. 
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Figure 10: Dahtz-plot density for the data. 
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I I F 




y y y 



Figure 11: Distribution over y. Points with errors - data, magenta histogram ~ signal {rrih — 
60MeV/c2,T,j = 10-i°s, \U^,h? = 10"^), red - K^2, blue - K^.^ and green - K^2'^. 



5.4- Simultaneous fit results 

As an example, the results of simultaneous fits in x-stripes 3 and 5 (m/i=60MeV/c^, 
Th = 10~^°s, Majorana type) are shown in Fig. [13], [HI 
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stripe 1 (0.2 <x< 0.25) Stripe 4 (0.35 < a: < 0.4) Stripe 7 (0.5 <x< 0.55) 




cosOn-t cosOa-y cosO, 



Figure 12: Distribution over cosO^^. Points with errors - data, magenta histogram - signal {rah = 
60MeV/c2,T,, = 10-iOs, |C/^/,|2 = lO'^), red - if^a, blue - and green - ^^27. 




Figure 13: Simuhaneous fit for m,,=60MeV/c2, t,, = lO^^^s. Stripe 3 (0.3 < x < 0.35). x^/n.d.f. = 
162.6/81. Points with errors - data, magenta - signal, red - KTr2i blue - if^a, green - K^2'y 




y cos6^^ 

Figure 14: Simultaneous fit for m;,=60MeV/c^ r/, = lO'^^s. Stripe 5 (0.4 < x < 0.45). x^/n-d.f. = 
188.3/81. Points with errors - data, magenta - signal, red - KTr2, blue ~ i^^a, green - K^2-f 

Both signal and background shapes are taken from MC. MC histograms are smoothed 
and the result is stored as f{z) function {z = y or cosOf^^). For better fit, we allow these 
functions to be slightly widen and shifted. We do it by using f{k-z + h) instead of 
f{z) in the fit, where fit parameters k and h are the same for signal and background and 
are different for y and cos9^^. For all selected x-stripes /c ~ 1 and 6 ~ 0, i.e. our MC 
describes the data properly (see [l^ for details). 

The simultaneous fit gives a signal event number in each x-stripe. As we use the same 
data several times we should take care of the correct estimation of a statistical error. The 
whole procedure of the simultaneous fit is as follows: 
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do simultaneous fit of two histograms and obtain {pi} - best parameter values (they 
correspond to the global minimum); 

take {pi} as initial values and perform x^/n.rf./. and error estimation for one his- 
togram cos 6^^ using a single call of MINOS program [16 . 



6. Signal efficiency 

A signal efficiency is defined as a number of events passed all cuts and fallen into the 
final kinematic region divided by an initial event number in the same kinematic region. 
Efficiency distributions for different lifetimes are shown in Figs. [15] -[17] as a function of 
Uh mass. 
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Figure 15: Signal efficiency e vs Figure 16: Signal efficiency e vs Figure 17: Signal efficiency e vs 



I'h mass for th = 10 ^s. 



Vh mass for Th = 10 s 



Vfi mass for th = 10 s. 



Two factors influence the signal efficiency: 

• Uh effective lifetime is larger for small rrih because of a Lorentz 7-factor; 

• a signal peak in y moves towards small muon energies for large rrih masses and is 
partly cut by y > 1 selection; it results in a lower efficiency for large Uh masses. 



7. Upper limit on 

7.1. \U^h\'^ extraction from X- stripes 

From the simultaneous fit we get Nexp{K — )■ fiUh) for each x-stripe. This event number 
should be transformed into \U^h\'^. This could be done either for each x-stripe or for the 
entire selected kinematic region. To avoid a systematics depending on an x-stripe (for 
example backgrounds could be described worse by MC in a certain x-stripe and shift a 
total result) we use the ffist method. 

As a normalization it is natural to use decay. The main formula is the following: 

Nexp{K^2'y) is taken from our previous analysis [l^ in a wide kinematic region, e{K — )■ 
^Vh) and e{Kn2j) are efficiencies obtained from MC. BR{K^2'y) is taken from the theory 
because an experimental measurement is very old and has a large error (the mean value 
is consistent with the theoretical prediction). In future, this BR could be measured using 
ISTRA+ data. 

BR{K — )■ fiUh) is substituted by the following expression: 
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Here BR{K^2) is taken from PDG [18|], f{mh) contains chirality flip and phase space 
factors. For the Dirac case from the general formula in we get: 

"^ft 1- ^ + 2^ + — ^ 1 f 2 2 V 

fD{mh) = — ■ \ i^ + zzT- zj-) - 4— f • 



m 



m 



K 




foiji^h) is dominated by the chirality flip factor 1 + ij^Y (see Fig. [T8|) and for 
the mass interval irih = 30-80MeV/c^ it varies from 1.1 to 1.6. For the Majorana case, 
fM{mh) = 2 ■ fD{mh). 

Finally for lU^hl"^ we get 
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Figure 19: Final fit, m/j^SOMeV/c^, t,, = IQ^^^a. Figure 20: Final fit, m/i^SOMeV/c^, r,, = lO^^^s. 



Dirac case. 



Majorana case. 
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1.2. Averaging |f/^/iP and setting upper limits 

Values |f/^/i,p are calculated for all seven x-stripes (we will denote them as 
and then averaged {\U^'j[\'^). We call this averaging procedure a final fit. The final fit 
for a certain {nih, Th) is shown in Fig. [19] (Dirac case) and Fig. [20] (Majorana case). The 
effect in (only a statistical error is considered here) does not exceed (1 — 2)cr for 

all {rrih, r^) and hence an upper limit should be set. The averaged value is used for a 
one-sided upper limit (U.L.) calculation: 

Tav 1 2 



Here autot is a total error of If/u/^P measurement 



U.L.{95%C.L.) = \U^l\' + 1.64 ■ au,tot. 

fj.h I 



7.3. Systematic error 

Main sources of the systematics are: 

• fit systematics; 

• cut on X (number of x-stripes in the final fit); 



• cut on y in x-stripes; 

• x-stripe width; 

• bin size in y and cos9^^ histograms; 

• cut on z^tx- 

The largest contribution to the total systematic error comes from the fit systematics 
caused by non-ideal MC shapes of the signal and the backgrounds. To estimate this 
systematic error the following procedure is used: 



errors for |f/*^''P are scaled for each x-stripe proportionally to ^J jn.d.f .\ 
averaging is repeated with these new scaled errors; 

new averaged value |[/"^''''="''^'^|2 j^g^g larger error a scaled which is treated as a scaled 



^stat + fit- Here astat is a statistical error of If/^"'^ 



'stat -r "^systjif ""^^^^ ^ stat ^ ciiwi 

The systematics of a cut on x is estimated as follows: 

• averaging is done for different number of x-stripes in the fit (varying cut on x); 

• the dependence of on x-cut is fitted by a straight line; 

• the slope of this line multiplied by the x-stripe width is the estimation of the 
systematic error. 



Details of this procedure are described in [19 



The systematics of a cut on y is calculated in a similar way by varying y-cut value 
and fitting the dependence of on y-cut by a straight line. 

x-stripe width is changed {dx = 0.035, dx = 0.07) and the whole procedure (simul- 
taneous fits in x-stripes, final fit) is repeated for new dx. New values of \U^h^^^\'^ are 
compatible with old ones and hence no systematics is found here. 

The systematics caused by the bin size in y and cos 9^^ histograms is estimated in a 
similar way and the result is the same: no additional error is found. 

The cut on z^tx is varied within the resolution. The systematical error is calculated 
similar to that of the cut on x . 
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l.J^-. Upper limits 

Upper limits for different lifetimes as a function of rrih are shown in Figs. [211 ESI The 
limits are calculated for the following values of mh. 30, 40, 50, 60, 70 and 80 MeV/c^. 
The curve in the figures is the interpolation between these values. The upper limits 
could be compared with the region predicted in ^ (shown with a blue stripe). As an 
example, contribution of all errors to the final result is shown in Table [1] for the Dirac 
case, Th = 10~^°s. 
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Figure 21: Upper limit for |?7p/j,P vs (Dirac case). Black line - obtained upper limits, blue stripe 
prediction from Q. 
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Figure 22: Upper limit for |?7^/ip vs mh (Majorana case). Black line - obtained upper limits, blue stripe 
- prediction from [6]. 



Exact numbers are collected in Tables |2] and |3l 
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<^stat 


^syst,x 


^syst.y 


^ syst.z 


^syst,fit 


U.L. 


30 


0.6 


0.4 


0.001 


0.002 


0.1 


0.5 


1.6 


40 


0.7 


0.4 


0.01 


0.009 


0.04 


0.5 


1.8 


50 


0.7 


0.4 


0.01 


0.004 


0.06 


0.5 


1.8 


60 


0.8 


0.4 


0.005 


0.03 


0.4 


0.5 


2.0 


70 


0.9 


0.5 


0.0002 


0.2 


0.4 


0.5 


2.2 


80 


1.0 


0.5 


0.08 


0.3 


0.7 


0.6 


2.8 



Table 1: Statistical and systematic errors for i/h of the Dirac type, — 10 ^"s. Fit results and C/.L.'s 
are in 10^^ units. 



8. Conclusions 

We have performed a search for a heavy neutrino of the Dirac and Majorana type in 
K — )■ fiUhi^h ^l) decay assuming that Uh is a part of z/^ flavor eigenstate and decays 
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rrih, MeV/c^ 




IQ-'^s 


U.L.,Th = 10-"'s 


U.L.,Th = 


10-"s 


30 


7.9 




1.6 


1.4 




40 


7.0 




1.8 


1.8 




50 


5.4 




1.8 


2.0 




60 


4.6 




2.0 


2.6 




70 


4.3 




2.2 


3.2 




80 


4.6 




2.8 


3.3 




Table 2: Upper limits for Vh of the Dirac type. U.L.'s are in 10 


units. 


nih, MeV/c^ 


U.L.,Th^ 


lO^^s 


U.L.,rh = 10-^"s 


U.L.,TH^ 


10-"s 


30 


4.5 




1.0 


0.8 




40 


4.4 




0.9 


0.9 




50 


2.9 




1.0 


1.0 




60 


2.2 




1.0 


1.4 




70 


2.2 




1.3 


1.2 




80 


2.2 




1.7 


1.1 





Table 3: Upper limits for of the Majorana type. U.LJs are in 10 ^ units. 



radiatively into a massless neutrino and a photon and obtained upper limits at 95% CL. 
for the mixing matrix element {U^hl"^. 

The upper limit at 95% CL. in a mass region SOMeV/c^ < < 80MeV/c^ for 
10-"s <Th< lO-^s is U.L. ~ (1 ^ 5) ■ 10-5 (Majorana type of Uh) and U.L. ~ (2 ^ 8) ■ 
10"^ (Dirac type). The obtained values close the allowed region for |f/^hP suitable for 
LSND/KARMEN/MiniBooNE anomaly explanation proposed in 

Authors would like to thank S.N.Gninenko, D.S.Gorbunov, V.A.Rubakov, A.A.Saratov 
(INR RAS) and R.Shrock (YITP, Stony Brook) for numerous discussions. The work is 
supported by the Russian Fund for Basic Research (grants 10-02-00330-a and 11-02- 
00870-a). 
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